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Abstract 

We study the three-dimensional XY-spin glass as a model for the resistive behavior of 
granular superconductors containing a random distribution of it junctions, as in high-T c 
superconducting materials with d-wave symmetry. The tt junctions leads to quenched in 
circulating currents (chiralities) and to a chiral-glass state at low temperatures, even in 
the absence of an external magnetic field. Dynamical simulations in the phase represen- 
tation are used to determine the nonlinear current-voltage characteristics as a function 
of temperature. Based on dynamic scaling analysis, we find a phase-coherence transi- 
tion at finite temperature below which the linear resistivity should vanish and determine 
the corresponding critical exponents. The results suggest that the phase and chiralities 
may order simultaneously for decreasing temperatures into a superconducting chiral-glass 
state. 
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Granular superconductors, including the high T c materials, are often described by arrays of 
superconducting grains coupled together by Josephson junctions. In a system of conventional 
junctions (s-wave pairing) the phases of the superconducting order parameter of neighboring 
grains tend to be locked with zero phase shift and a phase- coherence transition occurs for 
decreasing temperature into a state with long-range phase coherence and vanishing linear resis- 
tivity |lp. On the other hand, there is increasing evidence for <i-wave pairing symmetry in high 
T c materials with the remarkable consequence of the appearance of "71"" junctions (negative 
coupling) characterized by a phase shift of n. This leads to frustration effects in granular 
samples, even in zero external magnetic field, since a closed loop containing an odd number 
of 7r junctions gives rise to an spontaneous circulating current (orbital magnetic moment) and 
can be a possible explanation of the paramagnetic Meissner effect P, 0, [|. There is a close 
connection between this system and the XY-spin glass model [|| where the two-component 
spins, s = (cos(#), sin(#)), correspond to superconducting grains and the random ferro or anti- 
ferromagnetic interactions to the Josephson couplings. A chiral order parameter can be defined 
in the XY-spin glass model measuring the direction of circulating current (vortex) in the loops 
for the superconducting array. Based on this analogy and numerical results for the XY-spin 
glass in three dimensions [|J, it has been argued that the equilibrium low-temperature state is 
a chiral glass but with no long-range phase coherence and therefore not a true superconductor. 
However, recent results for XY-spin glass suggest that a spin glass (phase coherence) transition 
is possible |j| in addition to the chiral glass transition. These results however rely on simula- 
tions for the ground state defect energy and spin-glass order parameter which are inaccessible 
experimentally. Experiments often measure transport properties and it would be of interest to 
study directly the current- volt age behavior from dynamical simulations and examine the issue 
of the finite-temperature resistive transition. 

To study the current- volt age characteristics we assume a resistively shunted Josephson- 
junction model for the current flow between grains represented by the phases 9i on a cubic 
lattice and simulate the nonequilibrium behavior using the Langevin equations || 

Co^ + ^-E = -Jo E sin^ - 9 3 - A^) + £ rfy (1) 

where R Q is a shunt resistance, C Q is the capacitance to the ground, J Q > is the Josephson 
coupling, rjij represents uncorrelated thermal noise with < %(£) 2 >= 2k B T/R Q to ensure 
thermal equilibrium, and = or tc, represents the phase shift across the junction, with 
equal probability of and tc corresponding to the standard XY spin glass [|, |5], |^, 0] . We use 
units where K/2e — 1, R Q — 1, I a — 1 and set J R 2 C = 0.5, corresponding to the overdamped 
regime. The equations were integrated numerically || for a system of linear size L and the 
voltage V (electric field E/L) computed as a function of the driving current / (current density 
J = I/L 2 ) averaging over 5 — 10 different realizations of the A^ distribution. 

The nonlinear resistivity E/ J as a function of temperature T is shown in Fig. la, for L — 12. 
The behavior is consistent with a linear resistivity pl = limj^o E / J which is finite above a 
critical temperature T c . At lower temperatures, it appear to extrapolate to zero, indicating a 
superconducting transition with T c in the range 0.3 — 0.5. This is confirmed by a scaling anal- 
ysis of the nonlinear resistivity which assumes the existence of a continuous superconducting 
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transition at a finite temperature Jl] • Near the transition, measurable quantities scale with the 
diverging correlation length £ oc \T — T c \~ u and relaxation time r oc f , where v and z are 
the correlation-length and dynamical critical exponents, respectively. The nonlinear resistivity 
should then satisfy the scaling form |]] TE^ z ~ d+2 /J = g±{J£, d ~ 1 /T) where d is the system 
dimension and g(x) is a scaling function. The + and — signs correspond to T > T c and T < T c , 
respectively. A scaling plot according to this equation can be used to verify the scaling argu- 
ments and the assumption of an equilibrium transition at J = and also provide an estimate of 
the critical temperature and critical exponents. This is shown in Fig. lb, obtained by adjusting 
the unknown parameters so that the best data collapse is obtained. We estimate T c = 0.41(3), 
z = 4.6(4) and v = 1.2(4). The value of v agrees with the result obtained from calculations 
oi pi using the fluctuation-dissipation relation (zero current bias) in the vortex representation 
0. However, z is significantly larger but this could be a result of different dynamics. We 
also note that T c and the critical exponents agree with the estimates for the gauge-glass model 
H suggesting a common universality glass, although in the latter case Ay has a continuous 
distribution in the interval [0, 2ir]. Similar agreement has been found in the vortex represen- 
tation [0]. Since, in the phase model, the superconducting transition should correspond to the 
critical temperature where long-range phase coherence sets in ||, our results suggest a phase 
coherence transition at finite temperatures. This is in contrast with other numerical results 
indicating a phase coherence transition only at zero temperature However, it is consistent 
with more recent calculations showing evidence of a lower-critical dimension below or equal 3 
for the XY-spin glass model || and therefore a phase-coherence transition at finite temperature 
is possible. On the other hand, evidence of a chiral-glass transition has been found at T c ~ 0.32 
from Monte Carlo simulations ||], below our estimate for the phase-coherence transition, which 
could suggest an intermediate phase. However, this requires phase coherence in presence of 
uncorrelated chiral (vortex) disorder which is unlikely. Thus, the finite-temperature transition 
found in the present calculations should correspond to a single transition in the system where 
chirality and phase coherence order simultaneously into a low temperature superconducting 
chiral glass phase. Although our present conclusion is based on the scaling analysis of the 
nonlinear current- voltage characteristics which is a nonequilibrium property of the system, the 
good agreement with the previous results in the vortex representation obtained from calcula- 
tions of the (equilibrium) linear resistivity 0, suggests that this transition corresponds to the 
underlying equilibrium behavior. Additional calculations for the linear resistivity will help to 
settle this interesting issue. 
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Figure 1: Fig. 1. (a) Nonlinear resistivity E/J as a function of temperature, (b) Scaling plot 
near the transition obtained by adjusting T c , z and v where £ oc \T/T C — 
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